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AN  EXPERIMENT  IN  USING  THE  ELECTRONIC 
COMPUTER  "STRELA”  FOR  PROCESSING 
MEASUREMENTS  OF  METEORIC  PHOTOGRAPHS 

A.  K.  Sosnova 

The  processing  of  measurements  of  meteor  photographs  Is  very 
laborious . 

We  made  an  attempt  to  use  the  electronic  computer  "Strela" 
for  the  calculation. 

"Strela"  Is  directed  by  programming  and  automatically  solves 
problems  reduced  to  a  definite  sequence  of  arithmetical  and  logical 
operations*.  In  the  machine  numbers  are  represented  and  calculations 
with  them  are  carried  out  with  55  binary  symbols  ( significant  numbers) 
with  a  floating  decimal.  The  range  of  nvunbers  used  Is  +  . 

A  three-address  system  of  Instructions  and  a  parallel  means  of 
storing,  transmitting,  and  converting  numbers  and  Instructions  charac¬ 
terize  the  machine;  calculations  are  done  In  the  binary  system;  the 
conversion  of  the  Initial  numerical  data  from  the  decimal  to  the  binary 
system  and  the  reverse  conversion  of  the  output  Is  aoooflq)ll8hed 

*A.  I.  Kitov.  Electronic  Computers.  Isd.  "Sovetskoe  Radio." 

Moscow,  1936 • 
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automatically  by  the  machine. 

Before  proceeding  to  solve  any  problem  on  the  machine,  one  sets 
up  a  schedule  according  to  the  chosen  method  of  numerical  solution. 

The  numerical  method  Is  chosen  so  as  to  use  the  simplest  and  most 
uniform  formulas,  although  in  doing  ao  It  Is  often  necessary  to 
sacrifice  the  required  accuracy  for  an  Increase  in  the  scope  of  the 
calculations. 

The  method  must  reduce  the  solution  to  elementary  operations 
and  provide  control  of  the  solution. 

Into  the  flow  pattern  of  each  problem  there  enters  a  sequence 
of  arithmetical  operations,  of  provisory  and  absolute  transfers  of 
control,  of  change  and  renewal  of  instructions. 

The  complete  program  of  the  solution  is  formed  In  conformity  with 
the  flow  pattern  and  the  earlier  made  distribution  of  the  data  In  the 
registers  of  the  machine. 

The  program  for  solving  the  problem  written  In  the  form  of  a 
system  of  Instructions  and  the  Initial  numerical  data  are  entered  on 
punch  cards.  Each  number  or  Instruction  Is  recorded  on  one  line  of  a 
standard  80-column  punch  card  with  a  certain  combination  of  holes. 

The  numbers  and  Instructions  are  put  Into  the  codes  with  which 
the  machine  operates.  Each  code  contains  42  binary  digits. 

The  numbers  are  presented  In  semllogarlthmlc  form:  each  nximber 
N  Is  represented  in  the  form  of  2  groups  of  digits— a  mantissa  and  an 
ordinate 

Nio  -  d-loP,  0.1.^d<l, 

where  nsintlssa  ^  Is  the  sequence  of  the  digits  representing  the  nundber 
and  the  ordinate  2  Indicates  between  which  numbers  of  the  mantissa  the 
decimal  belongs. 
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Every  decimal  digit  of  a  number  N  In  the  machine  Is  represented 
In  the  binary  system  ( binary-decimal  numbers)  .  The  signs  of  the 
numbers  In  the  machine  are  plus-zero  and  minus-unity  in  the  nmber  set 
aside  for  the  sign. 

If  we  examine  the  code  as  an  Instruction,  It  consists  of  the 
following  parts:  1st,  2nd,  3rd  address,  control  signal,  and  number  of 
operation  to  be  carried  out. 

In  the  1st  and  2nd  addresses  the  Instructions  Indicate  the 
number  of  compartments  from  which  the  nvunbers  are  taken  on  which  are 
to  be  performed  the  operation  Indicated  by  its  own  number;  the  3rd 
address  Indicates  the  compartment  Into  which  the  number  obtained  from 
the  operation  Is  to  be  placed. 

Before  solving  a  problem  a  pack  of  punch  cards  with  the  program 
and  Initial  data  are  fed  Into  the  machine. 

The  work  of  the  machine  consists  In  sequential  execution  of  the 
Instructions  In  the  program. 

The  results  of  the  solution  are  fed  out  of  the  machine  onto  punch 
cards  and  are  printed  In  the  form  of  tables  of  numbers  In  the  order  pre¬ 
determined  by  the  program. 

In  the  FWU  "Strela  4"  electTOnlc  computation  center  we  composed 
and  debugged  a  program  for  processing  photographic  observations  of 
meteors . 

We  present  below  the  method  of  deriving  the  Initial  data  and 
the  flow  pattern  of  the  calculations  according  to  which  the  program 
was  composed. 

Derivation  of  Initial  data.  We  will  designate  the  point  of 
observation  at  which  the  shutter  was  fixed  as  A,  and,  correspondingly, 

we  will  call  the  negative  Interruptedly  tracking  the  meteor  negative 
A. 
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We  will  designate  a  corresponding  point  as  B  and  the  negative 
obtained  there  we  will  call  negative  B.  We  will  adopt  point  A  as  the 
base  point  and  In  respect  to  It  we  will  determine  the  heights  and 
velocities  of  the  points  of  the  meteor,  as  well  as  the  distance 
from  the  zenith,  and  the  radian  distance,  Zj^  (we  are  assvimlng  that 
the  negatives  are  made  In  guided  cameras)  . 

The  choice  of  reference  stars.  Reference  stars  are  selected, 
up  to  12  and  not  less  than  6,  uniformly  along  the  course  of  the  meteor 
and  on  bo  h  sides  of  it  on  each  negative,  and  one  star  not  farther  than 
1“  from  Its  optical  center.  (For  NAPA  3C/25  factory-produced  cameras 
we  may  take  the  Intersection  of  the  two  diagonals  of  the  frame  as  the 
optical  center  .) 

It  is  desirable  that  the  reference  stars  be  of  the  same  magnitude. 
It  Is  better  to  take  those  of  the  8th-9th  magnitude  as  references  for 
guided  cameras  and  those  of  the  6th-8th  magnitude  for  non-gulded  cameras. 

The  equatorial  coordinates  of  the  stars  are  written  out  for  the 
equinoctial  time  of  observing  the  meteor  fr(»n  AQKa  or  the  Yale  Catalogue 
with  an  accuracy  of  0.01"  for  a  and  O.l"  for  6  (the  time  Is  accurate 
to  0.1  of  a  year)  . 

Measurements.  The  plate  (or  the  film  frame  between  two  pieces 
of  glass)  Is  oriented  In  a  measuring  machine  so  that  one  of  the  axes 
of  the  device,  we  shall  call  It  the  x  axis.  Is  directed  along  the 
meteor's  trail. 

On  plate  A  are  measuz*ed  x  and  ^  for  the  reference  stars  and  the 
star  taken  as  the  optical  center,  for  l6  points  \uilformly  spaced  along 
the  trail  of  the  meteor,  and  for  the  ends  (or  mld-polnts)  of  the  gaps 
In  the  trail.  The  x  and  jr  coordinates  of  the  l6  points  are  subsequently 
used  for  finding  the  equatorial  coordinates  of  the  pole  of  the  meteoric 
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circle;  and  the  x  and  £  coordinates  of  the  ends  (or  mid-points)  of  the 
gaps,  for  determining  heights  and  velocities. 

On  plate  B  are  measured  the  x  and  ^  coordinates  of  the 
reference  stars,  of  the  star  taken  as  the  optical  center,  and  of  the 
l6  points  along  the  meteor's  track. 

In  order  to  avoid  the  effect  of  a  possible  movement  of  the  plate 
in  measuring  it  is  best  to  proceed  as  follows:  measure  first  the 
coordinates  of  the  stars,  then  of  the  points  along  the  meteor's  trail, 
aid  finally  the  coordinates  of  several  stars  again.  If  the  first  and 
last  measured  coordinates  of  the  same  stars  coincide  (within  the  limits 
of  measurement  error,  l.e.,  2-?m>)  we  may  consider  that  the  plate  did 
not  move. 

In  order  to  exclude  systematic  errors  the  measurements  are  made 
with  a  reversing  prism  or  by  turning  the  plate  180“. 

Two  applications  each  of  the  cross  hairs  of  the  measuring 
Instrument  are  made  for  each  measured  object  in  one  position  of  the 
prism  or  plate. 

For  our  calculationo  it  is  necessary  to  have  the  following 
values: 


a)  for  the  12  reference  stars  —  equatorial  coordinates  a^,  6^^ 
and  the  measured  coordinates  x^^,  y^^; 

b)  for  the  optical  center  stars  —  equatorial  coordinates  o^,  6^; 
the  measxired  coordinates  x^,  y^; 

c)  for  the  l6  points  along  the  meteor's  trail  —  the  measured 
coordinates  3^,  y^^ ; 

( values  a,  b,  and  c  are  necessary  for  both  plates  A  and  B) ; 

d)  for  the  20  points  at  the  ends  (or  mld-polnts)  of  the  gaps  In 

the  meteor's  track  on  plate  A  —  the  measured  coordinates  y^j^  ,  and 

the  times  t  cor3?e spending  to  each  of  these  20  points  (t  prefigured 
n  n 
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from  the  known  closing  speed  of  the  shutter) ; 

e)  there  must  also  be  given  the  latitude  qp  of  point  A  ( sin 

cos  qp)  and  data  on  the  base:  hour  angle  and  declination  t^,  of  the 
bearing  of  the  base  line  and  the  length  of  the  base  b  In  km;  and  for 
directional  cameras  the  sidereal  time  S  of  the  appearance  of  the  meteor 
(exactly  determined)  . 

For  a  given  base  of  observation  q),  b,  t^^,  and  6^  are  constant. 

The  program  of  calculations  Is  figured  for  the  12  reference 
stars j  l6  base  stars,  and  20  end  points  of  the  gaps  In  the  meteor's 
track. 

The  calculations  may  be  carried  out  also  for  a  smaller  niimber  of 
reference  stars  and  points  along  the  meteor's  trail  (but  this  number 
must  not  fall  below  6)  . 

In  addition,  five  more  values  must  be  given: 

f)  a  nvunber  one  less  than  the  number  of  reference  stars  for 

plate  A,  and  the  same  for  plate  B,  l^-l  (the  center-point  star 

Is  left  out  of  consideration) ; 

a  number  one  less  than  the  number  of  base  points  along  the 
meteor's  trail  on  plate  A,  k^-1;  and  the  same  for  plate  B,  k^-l; 

a  number  one  less  than  the  number  of  end  (mid)  points  In  the 
trail  gaps  on  plate  A,  m^-1. 

All  the  Intlal  data  enumerated  above,  beginning  with  Item  (a) 
are  put  on  punch  cards  and  these  are  fed  Into  the  machine  to  execute 
the  calculations. 

The  coordinates  of  the  poles  of  the  meteor  great  circles  A  and  B 
and  the  coordinates  of  the  radian;  the  heights,  velocities,  decelera¬ 
tions,  and  cos  (see  above)  for  the  20  points  are  calculated  from 
the  formulas  below. 
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Calcvilatlon  foxTHUlast 

1.  There  are  made  the  equalities  -  xo  ■  x^> 

.  ... .  ^ 


-  yo 
yv  -  yo 


ym“y° 


yi> 


where  i  ■  1  ...12,  k  «  1  ...l6,  m  -  1  ...20. 

2.  The  ideal  coordinates  of  the  reference  stars  are  figured: 


I, 


cos  sin  («^  — s„) 

'I  *  sin’t ”sln\+ cos  i'  cm  J.'itos  lo^ 

sln>j  cost,— cosS^  slnt,cos(a^  — «„) 
"siT»j"slnr»,+cos»^  co»:,  cos  (Sj  — s,)' 


(1) 


5.  By  the  method  of  least  squares  are  figured  the  coefficients 
of  the  relationship  between  the  measured  and  Ideal  coordiantes  for  the 
reference  stare  a,  b,  c,  d,  e,  f,  p,  s,  q,  t: 


—aXf  +byi  +c+pxi  *+sXi  - 

^  +ey,  +f+qx,  *+tXi  - 


(2) 


(3) 


For  a  check  the  machine  delivers  the  residual  errors  in  Eq.  2 
and  from  the  size  of  which  we  may  Judge  the  accuracy  of  the 

measurements  and  the  presence  or  absence  of  errors  in  identifying  the 
reference  stars. 

Prom  Eq.  2  and  are  found  by  substituting  Xj^  and  and 
from  them  Oj^  and  by  using  the  formulas: 

d/^  '-arcten^^ 

h 

(“i  *o)  —  •• ' 

co%{a^  — 

Then  from  and  by  the  method  of  least  squares  are  calculated 
the  coordinates  of  the  pole  of  the  meteor's  great  circle  A  and  sin  Rj^, 
where  is  the  declinations  of  the  points  and  from  the  calculated 
great  circle: 

cos  rtin  (^) 

iit-1  .  .  .  16, 
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sin  Rh  =sin  1),  sin  2^^+cos  cos  cos(a;i  ~*pjd‘ 


(4a) 


Similar  calculations  are  also  made  for  plate  B,  after  which  are 
determined  the  radian  coordinates  and  6^: 

tM  *4  Utapy^  >rimApA  sec  secaj;=-l;  |  ^  5^ 

tils  “i?  =— 1*  I 

The  following  Is  also  calculated 

cos  sin  t-cos  cos  2p^  cos  (^p^-~^pb), 

Where  Q  Is  the  angle  between  the  meteor  circles  A  and  B  at  the  radian. 

Finding  the  coordinates  for  simultaneous  points.  By  substituting 
Xm>  yj^  for  each  m  of  Eq.  2  are  calculated  and  from  them  a^, 

• 

Am 

Prom  Formulas  5a  —  the  analogue  of  Formulas  5  —  are  found  the 

auxlllaiy  values  tan  a  ,  tan  sec  a  : 

m  m  m 


tsBa„tsn«4  +tMiS^  secB^tiftt  secs*  =— 1; 

sec  sec  =  -1, 


(5a) 


where  =  S-t^j  from  5b  are  figured  the  coordinates  o^,  of  the 
point  simultaneous  to  6^: 


Sw  »m  *"•  *SP^Sm  secajn,=  -l;  | 

tsB  secips^B^  sec  I 


(5b) 


The  distance  of  point  m  from  point  A  Is  determined  from  the 
formula 


icosSj  sin  ) 


and  the  height  of  this  point  relative  to  the  horizontal  plane  of 
point  A 

"m  obs.  ■  'm'  t®’ 
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where  co»  =.siii  stn  ?--co8 cos  ?  cos  (5— */)„). 


(7) 


The  zenith  dletance  of  the  radian  Is  also  figured. 

Thus  by  combining  the  Initial  data  and  the  results  of  our 
calculations  we  have  a  number  of  values  for  and 

corresponding  to  them. 

Velocities  and  decelerations.  The  height  presents  Itself 

In  the  form  of  a  polynomial  of  the  4th  degree: 

Hobg  =■  a  +  bT  +  CT*+  pT®+ST*.  (9) 

The  parameters  a,  b,  c,  p,  s  for  ^^jg  (m  »  1...20)  are 
discovered  by  the  method  of  least  squares. 

Prom  a,  b,  c,  p,  s  are  calculated  the  values  of  calc 
subtractions  ^^g  -  calc  made;  for  all  values  of  nj  the 

velocities  and  decelerations  are  figured  by  differentiating  Formula 
9  with  respect  to  t: 

-r\2sx*„). 

The  machine  delivers  cos  z,^,  , 

«m  obs.^  \  calc.'  V  V  V  residual  declinations  of  Eq.  2 

cos  the  values  of  a,  b,  c,  p,  s  from  Eq.  9 
divided  by  cos  z^;  Blna^*,  sin  R^. 

As  a  check  on  the  calculations  they  are  carried  out  twice.  The 
approximate  amount  of  time  spent  on  processing  the  photographs  of  one 
meteor  used  as  a  basis  using  the  "Strela”  ccnnputer  Is  as  followst 
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choice  of  stars  and  transcribing  their  positions  from 


the  catalogues  4  hrs 

measurement  2  hrs 

preparation  of  the  Intlal  data  for  pvinchlng 
(filling  the  cards)  1  hr  30  min 

punching  and  checking  the  cards  30  min 

calculation  2.5  min 

printing  the  results  3  min 

Total  8  hr  6  min 


Thus,  preparing  the  Initial  data  for  the  calculations  and  the 
calculations  themselves  on  "Strela”  take  a  little  more  than  2  hours. 

At  the  same  time  no  less  than  50  hours  are  spent  on  processing 
the  measurements  when  calculating  by  hand. 

Using  a  "Strela"  computer  In  the  spring  of  1958  part  of  the 
photographic  meteor  observations  of  the  Stallnabad,  Kiev,  and  Odessa 
observatories  were  processed  in  the  MOU  Computation  Center. 

I  consider  it  my  pleasant  duty  to  thank  Prof.  Ye.  Ya.  Bugoslavskaya 
and  Prof.  G.  S.  Roslyakov  whose  advice  helped  me  very  much  In  carrying 
out  the  present  work. 

Supplement  1.  The  flow  pattern  and  the  progreunmlng  of  the  calcu¬ 
lations  are  set  up  In  such  a  way  that  they  permit  us  also  to  process 
negatives  gotten  from  non-gulded  ceuneras  —  It  Is  only  necessary  to  Im- 
sert  the  following  on  the  Initial  data  sheets:  zero  on  the  S  line, 

-t|j  on  the  t^  line,  and  the  hour  angles  to*  ti...txa  in  place  of  the 
right  ascensions  Oq,  ai...aia  of  the  reference  stars. 

Supplement  2.  Flow  pattern  for  transferring  nxaibera  to  the  base 
60  to  those  to  the  base  2  (60-*  2' 

"Strela"  conducts  trigonometric  operations  on  angles  expressed 
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abstractly  and  to  the  base  2.  Since  the  initial  data  a  and  6  are  given 
in  the  sexagesimal  system  —  in  hours,  minutes,  seconds  of  time  and 
degrees,  minutes,  seconds  of  arc  —  it  was  necessary  to  find  some  means 
of  changing  a  and  6  into  abstract  measure  in  the  binary  system.  The 
following  algorithm  was  used  in  our  program  for  that  purpose. 

If  the  angles  are  written  in  a  semllogarlthmlc  form  placing  the 
ordinate  02  before  all  numbers  (angles) ,  then  mantissa  d  of  angle  a, 
expressed  in  degrees  and  parts  of  degrees  —  minutes  and  seconds  of  arc 
(or  hours,  minutes,  and  seconds  of  time)  -  may  be  expressed  as 

=  10"‘  +/>j  />, 

+  ur“»r-f  A. 

where  bi,  b2...are  sexagesimal  divisions  of  the  nianber  (angle),  counting 
from  left  to  right:  bi  represents  tenths  of  degrees,  b2  —  units  of 
degrees,  ba  —  tenths  of  minutes,  etc. 

In  confomlty  with  this,  for  angles  expressed  both  in  arc  and  time, 
when  transposing  6o-*2  we  used  the  algorithm 

(  .  . .  ( X 1  lO+Xa)  lO+Xa)  6+X4)  IO+X5)  6+Xe)  IG+X7)  10+Xb) 

104x9)  10“*6”*2^,  where  «  ^1^~*  ^1  1-th  decimal  figure  of 

the  transposed  number  (the  angle  is  written  in  the  form  of  a  decimal 
number.')  . 

For  transposing  angles  into  abstract  measure  they  are,  as  usual, 
expressed  in  degrees  and  multiplied  by  the  multiplier  ^  and  when 
expressed  in  hours,  ^  ^  • 

The  60-»2  transposal  program,  included  in  the  complete  computa¬ 
tional  program,  enables  us  to  introduce  the  initial  data  for  the  angles 
directly  into  the  machine  in  hours  ( degrees) ,  minutes,  seconds,  and 
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fractions  of  seconds,  written  as  follows: 


+0.234852425  •  10= 

4-0,000001208  •  10* 
+0.6538171  •  10» 

-(VXI87I34  •  10* 


Instead  of  ,  03*  4b'"  52*  ,  425 
00*  00'"  01*  .208 

+65‘’38'17'’,1 

-00*2n3*,4 
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DisrarsuTioN  lisj 


OF  D3F3USB 


I'&r.  Copies  AIR  COHLIAIOS  Nr.  Copio 


i;i:;oQ,u.mi3HS  usap 

APCI:J-3R2  1 

iulL  (ARR)  1 


CTRRR  AGENCIES 


0  i  A  1 

i>  iwi*  ^ 

AlA  •  2 

Oi'S  2 

ACC  2 

P>.3  1 

NASA  i 

HAI®  1 


AP30 

SCPTR  1 

ASTIA  25 

ro-Bla  5 

■m-Blb  3 

ABDC  (ABY)  1 

SSD  (SSP)  2 

APGC  (PGP)  1 

BSD  (BSY)  1 

HADC  (RAY)  1 

APSWC  (SWP)  1 

APMTC  (MTW)  1 
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